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NOTATION 

1  bo 

Stratification  parameter  =  —  r** 

p0iy 

Rate  of  wake  growth  vertically  at  x  ■  0 

Initial  constant  rate  of  wake  growth  in  the 
vertical  direction  at  x  *  0 

Constant  in  empirical  formulae 

Observed  width  of  circular  Jet 

Gravitational  acceleration 

Constant  in  empirical  formulae 

Fluid  density 

Vertical  density  gradient 

Fluid  density  at  elevation  of  center  of  wake  or 
that  of  mixed  fluid  inside  the  wake 

Time  after  agitation  of  pendulum 

Time  of  collapse 

Time, after  time  of  collapse 
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in  pure  water 

Turbulence  intensity 
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tj  Mean  velocity  of  advance  . 

m 

U  Maximum  velocity  at  axis  of  Jet 

IT18a 

x,y  Coordinate  axes;  x-axis  is  horizontal  and  y-axis 

is  vertically  upward  with  origin  at  center  of 
wake 

xq  Horizontal  width  at  time  of  collapse 

yQ  Vertical  extent  of  wake  boundary  at  t  «  0 


y 

max 


Maximum  vertical  extent  of  wake  boundary  at  x  =  0 

and  t  =  t  - 
col 
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ABSTRACT 

V  \ 

The  force  of  gravity  causes  a  turbulent  wake  in  a  density- 
stratified  fluid  to  eventually  cease  its  vertical  growth  and  then 
to  collapse  towards  its  horizontal  midplane.  In  the  present  in- 
vestigation  this  phenomenon  has  been  Studied  experimentally.  The 
turbulent  wake  was  created  by  means  of  a  spiral  paddle,  agitated 
by  a  pendulum- type  arrangement  outside  a  transparent  lucite  tank. 
Data  were  obtained  from  tracings  of  the  motion  pictures  taken  by 
a  16  mm  movie  camera.  Both  the  pendulum  arrangement  and  the 
paddle  diameter  were  varied  to  find  the  possible  influence  of  the 
experimental  conditions. 


It  was  observed  that  the  initial  rate  of  growth  in  the  verti¬ 
cal  direction  is  constant,  depending  primarily  on  the  density 
gradient  and  the  agitation  mechanism  (i.e.  pendulum  and  paddle 
diameter).  This  initial  rate  of  growth  of  the  wake,  the  maximum 
vertical  thickness  of  the  wake,  the  time  at  which  collapse  begins 


and  the  turbulence  intensity  within  the  wake  at  that  time,  were 
all  correlated  with  the  Vaisala  frequency,  resulting  in  three 
important  constants  which  seemed  to  be  independent  of  the  experi¬ 


mental  conditions./ 

During  the  collapse  the  same  three  stages  of  collapse  could 
be  distinguished  as  had  been  defined  previously  by  Wuj  some  dif¬ 


ferences  in  the  post-collapse  behavior  were  observed,  however. 
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INTRODUCTION 

The  wake  of  a  body  moving  in  a  fluid  with  a  vertical  density 
distribution  is  considerably  different  than  that  of  the  same  body 
moving  in  a  fluid  having  no  density  gradient*  In  the  case  of  a 
density  gradient,  the  initial  expansion  of  the  mixed  fluid  is 
quickly  followed  by  a  collapse  in  the  vertical  direction,  which 
is  accompanied  by  a  further  spreading  in  the  horizontal  direction. 
This  phenomenon  is  caused  by  the  force  of  gravity.  The  volume  of 
the  fluid  in  the  wake  behind  the  body  has  a  more  or  less  constant 
density  due  to  mixing;  driven  by  the  hydrostatic  pressures  out¬ 
side  the  wake,  the  constant  density  fluid  is  forced  to  seek  its 
own  density  level  in  the  surrounding  fluid. 

Only  recently,  several  investigators  have  carried  out  labora¬ 
tory  Investigations  in  order  to  understand  the  interaction  of  the 
turbulent  three-dimensional  wake  with  a  linear  vertical-density 
gradient.  Schooley  and  Stewart  (lieference  l)  first  studied  this 
phenomenon  in  1962  with  a  model  self-propelled  body  and  measured 
the  vertical  collapse  and  the  horizontal  spreading  by  introducing 
dye  into  the  wake  behind  the  body.  Wu  (Reference  2)  in  1965 
studied  certain  aspects  of  the  wake  in  density-stratified  fluid 
after  collapse.  He  considered  only  gravitational  and  inertial 
effects  and  obtained  in  that  way  a  good  basic  understanding  of  the 
collapse  and  spreading  of  the  wake.  Moreover  he  analyzed  the 
internal  waves  generated  by  the  wake  collapse  phenomenon  (Ref¬ 
erence  6). 
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In  1966,  Stockhausen,  Clark,  and  Kennedy  (Reference  3)  in¬ 
vestigated  the  momentum! ess  three-dimensional  wake  in  a  stratified 
fluid  with  a  linear  density  gradient.  Their  experiments  were 
also  conducted  with  a  self-propelled  model  in  a  large  tank.  Varia¬ 
tions  in  salinity  were  measured  in  the  zone  of  mixing  by  means  of 
specially  designed  conductivity  probes.  In  this  way,  normalized 
isochlor  distributions  for  transverse  cross-sections  at  various 
distances  behind  the  model  were  deduced,  from  which  information 
could  be  obtained  about  the  .extent  of  the  mixed  region. 

In  the  present  study,  particular  attention  is  given  to  the 
growth  of  the  turbulent  wake  in  the  density-stratified  fluid  up 
to  the  moment  of  collapse.  For  some  of  the  runs,  information  is 
obtained  after  collapse  in  order  to  check  the  validity  of  some  of 
Wu's  findings  (Reference  2). 

EQUIPMENT  AND  EXPERIMENTAL  TECHNIQUE 
Tank  and  Stratified  Fluid 

The  experiments  were  carried  out  in  a  transparent  tank  48 
inches  deep,  43  inches  long,  and  12.5  inches  wide. 

The  vertical  density  stratification  was  obtained  by  intro¬ 
ducing  salt  water  through  a  diffusor  at  the  bottom  of  the  tank. 

A  plastic  screen  cloth  assured  a  uniform  outflow  from  the  diffusor. 
In  order  to  prevent  air  bubbles  from  disturbing  the  stratifica¬ 
tion  during  filling,  a  1/4-inch  brass  tube  was  mounted  at  the  top 
of  the  diffusor,  allowing  the  air  bubbles  to  pass  through  the  wall 
outside  the  lucite  tank  (See  Figure  l).  The  density  stratifica¬ 
tion  was  controlled  by  mixing  different  amounts  of  sodium  chloride 
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in  water.  The  successive  layers  of  salt  water  were  1  inch  thick 
and  alternately  colored  with  red  dye.  In  this  way,  24  layers  were 
brought  into  the  tank.  The  idealized,  initial  vertical  density 
distribution  of  the  fluid  in  the  tank  is  a  stepped  line;  but,  as 
has  been  verified  by  Wu  (Reference  2),  leaving  the  fluid  over¬ 
night  assures  a  linear  density  stratification  through  molecular 
diffusion. 

Paddle  Mixer 

The  turbulent  wake  was  generated  by  means  of  a  spiral  paddle 
which  was  known  to  generate  good  small  scale  turbulence  (see 
Reference  4).  A  spiral  paddle,  12-3/8  inches  long,  with  a  diam¬ 
eter  of  2  inches  was  first  used.  It  consisted  of  a  brass  spiral 
frame  (one  turn  in  1.5  inches)  on  which  a  screen  was  soldered. 

The  frame  was  fixed  in  place  by  three  discs  (see  Figure  2).  The 
paddle  was  supported  across  the  tank,  20  inches  above  the  bottom 
of  the  tank,  by  a  1/4-inch  O.D.  brass  tube,  in  such  a  way  that  the 
axis  of  rotation  of  the  paddle  is  perpendicular  to  the  front  and 
back  wall  of  the  lucite  tank.  In  the  brass  tube  many  small  open¬ 
ings  were  drilled.  A  rod  l/8-inch  in  diameter,  which  had  been 
painted  heavily  with  blue  dye  dissolved  in  polyvinyl  alcohol, 
was  inserted  into  the  brass  tube  immediately  preceding  a  test. 

It  took  several  minutes  for  the  dye  to  dissolve,  during  which 
time  the  turbulence  created  by  the  insertion  of  the  rod  was  al¬ 
ready  damped  out  as  could  be  clearly  observed  by  the  manner  in 
which  the  heavier  blue  dye  came  out  of  the  openings  of  the  brass 
tube.  The  paddle  was  then  agitated  by  means  of  a  pendulum-type 
arrangement  at  the  back  of  the  tank.  For  the  first  set  of  experi¬ 
ments,  this  pendulum  nad  an  effective  arm  of  21  inches  and  a  weight 
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of  1.5  pounds.  After  release,  the  pendulum  described  an  arc  of 
approximately  270°,  and  on  the  way  back  slightly  less  than  that. 
For  the  second  set  of  experiments  a  shorter  pendulum  was  used. 

This  pendulum  had  an  effective  arm  of  1  foot  and  a  weight  of  1.26 
pounds.  The  third  series  of  experiments  were  carried  out  with  the 
long  pendulum  and  a  spiral  paddle  having  a  diameter  of  4  inches 
(see  Figure  2). 

Data  Reduction 

The  wake  created  in  the  above-mentioned  way  was  more  or  less 
two-dimensional,  at  least  until  the  moment  of  collapse.  The 
boundary  of  this  wake  was  clearly  visible  against  the  background 
of  the  color-layer  stratification  system.  The  growth  and  collapse 
of  this  wake  was  photographed  by  a  16  mm  movie  camera.  With  the 
help  of  coordinate  lines  drawn  on  the  wall  of  the  tank,  wake  pro¬ 
files  could  be  traced.  The  time  could  be  read  accurately  to 
within  l/l 00  of  a  second  from  a  clock  which  was  also  photographed. 
A  set  of  sample  profiles  for  the  wake  generated  with  the  2-inch- 
diameter  paddle  and  the  long  pendulum  is  given  in  Figure  3» 

It  was  a  disadvantage  that  the  change  in  wake  thickness  after 
t  *  0  was  small  in  comparison  with  the  height  of  the  wake  at 
t  «  0.  Moreover,  irregularities  in  the  wake  shape  due  to  the 
presence  of  eddies  increased  the  difficulty  of  taking  accurate 
measurements,  and  resulted  in  a  scatter  of  the  data  points  along 
the  course  of  the  curves  (see  Figures  5>  6  and  7). 
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RESULTS 


Phase  I  of  this  study  was  performed  with  the  2-inch-diameter 

spiral  paddle,  agitated  by  the  long  pendulum.  Nine  tests  were 

carried  out  and  analyzed:  one  in  pure  water  and  eight  with  various 

linear  density  gradients.  The  particular  method  of  introducing 

the  salt  water  in  layers  allowed  a  minimum  value  for  a  ■  —  of 

PQ  dy 

0.002. 


In  pure  water,  where  the  turbulent  wake  is  more  or  less 
circular  in  cross-section,  the  area  of  the  wake  was  determined 
from  tracings  of  the  movie-film  by  means  of  a  planimeter.  Prom 
the  measured  area,  the  root-mean-square  of  the  radius  of  the  wake 
was  computed  and  plotted  versus  time,  as  in  Figure  5.  Time  t  -  0 
is  thereby  taken  as  the  instant  the  movement  of  the  pendulum 
stopped  completely, which  could  be  determined  from  the  film.  The 
curve  shows  a  continuously  decreasing  rate  of  wake  growth  in  pure 
water;  initially  the  rate  of  growth  Is  almost  constant. 

Tracings  were  also  drawn  from  the  boundary  of  the  wake  in 
the  tests  with  the  linear  density  gradients.  With  the  aid  of  these 
tracings,  the  vertical  ordinate,  y,  of  the  wake  boundary  at  x  *  0 
was  determined.  For  a  definition  sketch  of  the  coordinate  axes, 
see  Figure  4.  The  vertical  ordinates,  y,  were  plotted  versus  time 
in  Figure  5/  where  the  different  tests  are  identified  by  the  value 
for  a  ■*  ~  5?*  .  Although  there  occurs  a  certain  amount  of  scatter 

p0 

due  to  the  way  in  which  y  is  determined,  all  curves  show  the  same 
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tendency.  Not  only  does  the  rate  of  growth  in  the  vertical  direc¬ 
tion  diminish,  but  the  wake  height  reaches  a  maximum  and  decreases 
afterwards.  It  is  readily  3een  that  as  the  density  gradient  in¬ 
creases  (i.e.,  increasing  value  for  "a"),  the  rate  of  growth  of 
the  wake,  the  maximum  wake  height  y  and  the  time  t  _ ,  in 
which  this  maximum  height  was  reached,  all  decrease. 

In  order  to  find  the  possible  influence  of  the  pendulum  ar¬ 
rangement,  Phase  II  of  the  study  was  conducted.  The  arm  of  the 
pendulum  therefore  was  shortened  to  approximately  one-half  of  the 
length  used  in  Phase  I.  Consequently,  the  time  of  agitation  de¬ 
creases  by  about  10  percent.  With  this  new  arrangement,  five 
tests  were  performed:  one  in  pure  water  and  four  with  various 
linear  density  gradients.  The  results  of  Phase  II  are  given  ir. 
Figure  6. 

Since  the  paddle  size  could  also  have  an  effect  on  the 
phenomenon,  a  paddle  with  a  4-inch  diameter  was  built  and  six 
tests  were  carried  out  with  this  paddle,  using  the  long  pendulum 
(Phase  III).  Figure  7  shows  the  results  of  this  phase  of  the 
study.  Curves  of  the  same  type  are  obtained,  with  a  larger  ini¬ 
tial  height  of  the  wake,  y  . 

ANALYSIS  OF  DATA 

A  closer  look  at  Figures  5,  6  and  7  immediately  reveals 
certain  features  of  the  growth  of  the  turbulent  wake  in  the  ver¬ 
tical  direction.  It  is  apparent  from  these  curves  that  in  the 
majority  of  the  tests,  the  rate  of  growth  during  the  initial  stage 
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is  almost  a  constant,  indicated  by  the  'lope,  bi,  of  the  straight 
line  fitted  through  the  first-measured  points.  The  actual  values 
of  bi  for  all  tests  are  given  in  Table  1.  With  the  aid  of  bi,  the 
data  can  be  replotted  in  a  dimensionless  way,  which  gives  a  better 

y-y0 

picture  of  what  is  actually  taking  place.  When  plotting  — — 

yo 

versus  —^1.*  t  •  the  straight  line  sections  coincide  on  a  1:1 
*o 

slope  (see  Figures  5a>  6a  and  7a)*  Depending  on  the  density 

yo  r— 

gradient,  expressed  non-dimensionLlly  in  this  case  as  —  y  ag, 

the  individual  curves  branch  off  after  some  time,  showing  more 

clearly  that  the  maximum  height,  y  ,  increases  with  a  decreasing 

max 

density  gradient.  Moreover,  the  time  of  collapse,  t  ^  increases 

with  decreasing  density  gradient.  To  find  out  more  about  the  re¬ 
lationship  between  the  maximum  height  and  the  density  gradient. 


y  -  y 

'max  o 


is  plotted  versus 


bi 

y O  V~aS 


for  all  available  tests. 


The  result  is  given  in  Figure  8.  The  relationship  is  linear  and 
the  slope  is  equal  to  0.8.  This  means  that 


- -  =  1.25 

-  yo>  V«S 


[i] 
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for  all  tests,  independent  of  pendulum  length  or  weight,  paddle 
diameter,  or  density  gradient.  This  relation  is  especially  im¬ 
portant  because  now  it  is  possible  to  extrapolate  this  line  to 
values  of  "a"  which  occur  in  other  circumstances. 


ymax  "  yo 

Another  linear  relationship  can  be  found  if  -  is 


plotted  versus  t  —  (see  Figure  9).  The  slope  of  this  line  is 

coi  yQ 

0.64,  which  means  that 


fa  '  70  =  0.64  [2] 

*001 •  bl 

Combining  the  results  of  [l]  and  [2],  another  important  constant 
is  obtained: 


tool  =  1-25  [3] 

Equation  [3]  could  also  have  been  obtained  if  t  ^  •yj  ag  is  com¬ 
puted  for  all  tests  and  plotted.  Figure  10  gives  this  plot  and 
the  average  value  of  t  1  'yfag  is  again  1.25.  It  can  be  seen 

that  there  is  a  certain  amount  of  scatter,  but  it  is  difficult 
to  precisely  determine  t  ^  in  the  experiment. 
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DISCUSSION 

Growth  of  the  Wake 

In  order  to  obtain  a  better  understanding  of  what  is  actually 
taking  place  during  the  growth  period  of  a  turbulent  wake  in  a 
density -stratified  fluid,  the  generation  of  this  wake  in  the  pres¬ 
ent  experiments  will  be  discussed  first.  As  soon  as  the  pendulum 
starts  its  swinging  motion,  a  large  scale  eddy  is  created  which 
is  circular  in  shape  (see  Figure  11a).  This  eddy  grows  in  size 
through  the  diffusion  of  turbulence  within  it.  Although  the  shape 
of  the  wake  is  initially  almost  circular,  the  wake  flattens  as 
buoyancy  starts  playing  a  role.  In  fact,  at  t  *  0,  when  the  pen¬ 
dulum  stops  Its  motion,  the  wake  shape  is  not  quite  circular  (see 
Figure  lib),  perhaps  because  of  the  more  effective  turbulent  dif¬ 
fusion  in  the  horizontal  direction.  The  wake  still  grows  in  the 
vertical  direction,  but  slightly  less  rapidly  than  in  the  hori¬ 
zontal  direction.  The  energy-containing  eddies,  which  cause  the 
growth  of  the  wake  by  their  entrainment  of  surrounding  fluid,  are 
being  suppressed  more  and  more  near  the  top  and  bottom  of  the  wake. 
Finally,  at  the  time  of  collapse,  the  vertical  height  and  the  den¬ 
sity  difference  inside  and  outside  the  wake  become  so  large  that 
the  stabilizing  gravity  force  is  strong  enough  to  suppress  the 
turbulence  completely  in  the  vertical  direction.  At  this  moment 
there  is  only  growth  in  the  horizontal  direction  (see  Figure  11c). 
After  collapse,  the  gravity-induced  forces  become  predominant  and 
the  mixed  fluid  tends  tc  flow  sideways  to  the  plane  of  equilibrium. 
The  motion  of  the  wake  profile  thereafter  depends  on  the  degree 
of  mixing  and  the  density  gradient  (Figure  lid). 
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In  the  following  section,  the  instant  of  collapse  and  the 
horizontal  spreading  will  be  considered  in  detail. 

> 

Instant  of  Collapse 

A  rough  explanation  of  the  growth  process  of  the  wake  in  the 
vertical  direction  can  also  be  attempted  in  terms  of  turbulent 
energy.  A  water  particle  near  the  boundary  of  the  turbulent  wake 
moves  as  part  of  a  turbulent  eddy  or  puff.  Considering  only  its 
vertical  component,  its  kinetic  energy  per  unit  volume  can  be 
described  as 


K.E.  =  p  u'2  [4] 

o 

where  pQ  is  the  density  of  the  midplane,  assuming  that  the  fluid 

inside  the  wake  is  fully  mixed,  and  ~\J  u'2  is  the  turbulence 
intensity  in  the  vertical  direction. 

Because  of  the  presence  of  a  density  gradient  outside  the 
wake,  the  same  particle  becomes  surrounded  by  lighter  fluid  when 
crossing  the  boundary.  Assuming  again  a  completely  mixed  wake, 
the  density  difference  at  the  wake  edge  can  be  computed  as 


[5] 
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Wit  h  respect  to  the  midplane,  this  particle  thus  has  a  potential 
energy  expressed  as 


[6] 


At  the  instant  of  collapse,  when  the  motion  in  the  vertical  di¬ 
rection  has  stopped,  the  kinetic  energy  of  the  fluid  near  the 
wake  boundary  is  assumed  to  be  proportionate  to  the  potential 
energy  gained  due  to  displacement  from  the  midplane.  This  rejults 
in  the  following  expression: 


K.E. 

P.E. 


2  d£ 
max  by 


g 


constant 


[7] 


This  can  be  further  reduced  to 


constant 


[8] 


For  the  present  experiments,  ymax  has  been  measured  for  a 


certain  V~ag.  To  evaluate  the  turbulence  intensity  uq '2  . 
however,  a  little  more  has  to  be  said  about  the  turbulence  in¬ 
tensity  in  an  axially  symmetric  wake  in  pure  water. 


HYDRONAUTICS,  Incorporated 


-13- 

A  fully  turbulent  wake  is  separated  from  the  surrounding 
nonturbulent  fluid  by  an  irregular  turbulent  front  in  which  large 
scale  eddies  convect  the  fully  turbulent  fluid  outward  and  simul¬ 
taneously  bring  nonturbulent  fluid  inward.  Thus  the  large  scale 
eddies  control  the  rate  of  spreading  of  the  turbulent  wake,  while 
the  small  scale  eddies  bring  the  entrained  fluid  to  a  turbulent 
state.  Because  of  dissipation  and  energy  transfer  to  the  en¬ 
trained  nonturbulent  fluid,  the  turbulence  level  inside  the  wake 
will  lower  continuously.  Prom  the  principle  of  Reynolds  number 
similarity,  the  velocity  of  advance  of  the  turbulent  front  is 
proportional  to  the  root-mean-square  turbulent  velocity  : 


1 


constant 


[9] 


Townsend  (Reference  5)  has  verified  this  equation  and  found  the 
constant  to  be  O.51.  The  constant  will  be  approximately  the  same 
for  an  axially  symmetric  wake,  so 


dt 


[10] 
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Let  us  now  return  to  the  present  study  where  there  is  also 
the  effect  of  stratification.  In  order  to  find  a  suitable  value 

of  — —  at  the  instant  of  collapse,  let  us  try  to  understand 

the  behavior  of  the  mixed  region  in  the  turbulent  wake  during  the 
initial  stage  of  wake  growth.  In  the  early  wake,  turbulent  mixing 
is  the  dominant  feature  of  the  motion,  and  the  wake  expands  in  all 
directions.  If  a  particle  is  considered  at  some  x  and  y  (see 
Figure  4),  then  the  mean  turbulent  motion  is  almost  radial.  This 
radial  turbulent  vector  has  a  y-component  vertically  upward  for 
the  particle  shown  in  Figure  4  and  a  horizontal  component  in  the 
-x  direction.  Until  the  instant  of  collapse,  the  x-component  is 
affected  very  little  by  the  density  gradient,  while  the  upward 
motion  is  suppressed  by  the  gravity- induced  force  due  to  the 
density  gradient.  That  is,  a  particle  on  the  x-axis  is  not  sig¬ 
nificantly  influenced  by  gravitational  effects  during  the  period 
of  vertical  wake  height  increase.  During  the  collapse  of  the  wake, 
however,  the  gravity  forces  start  to  influence  the  behavior  of  the 
mixed  region  by  returning  the  fluid  containing  the  displaced 
salinity  toward  a  level  of  density  equilibrium  where  it  spreads 
laterally.  This  means  that  during  the  collapse  stage  a  change 
in  the  horizontal  spreading  is  to  be  expected.  Before  collapse, 
however,  the  horizontal  spreading  should  be  the  same  as  in  the 
pure  water  case.  This  concept  will  be  checked  now  with  the  aid 
of  experimental  results.  In  Figure  12, 
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x-x 


versus  time  for  the  pure  water  case  and 


versus  time  for 


several  runs  in  density-stratified  fluid  are  plotted.  The  plots 
show  that  indeed  in  the  beginning  there  doesn't  appear  to  be  much 
difference  between  the  wake  growth  in  fresh  water  and  that  in 
density-stratified  fluid.  At  the  moment  of  collapse,  the  line  is 
dotted  and  the  slope  of  the  line  is  very  close  to  the  slope  in 
the  pure  water  case.  The  curvature  changes  after  collapse,  as  can 
be  clearly  seen  in  Figure  12.  This  is  due  to  the  sudden  lateral 
surge  when  the  heavier  fluid  flows  down  and  sideways  to  its  own 
density  level.  Note  that  the  change  in  curvature  is  larger  and 


earlier  with  increasing  a 


-i- |£  as  was  to  be  expected.  Since 


it  has  now  been  shown  that  the  turbulence  intensity  that  would 

have  existed  in  the  vertical  direction  at  y  and  at  t  .  in  the 

J max  col 

case  of  density -stratified  flow  can  be  found  by  calculating  the 

slope  of  the  pure  water  v?  versus  time -curve  at  tCQl,  the  ratio, 
in  Equation  [8],  can  be  computed  for  each  run.  These  values  of 


ymax  V* 


are  tabulated  in  Table  1  and  appear  to  be  more  or  less 


constant,  with  an  average  value  of  0.36.  In  general,  then,  it  can 
be  stated  that  approximately 


HYDRONAUTICS,  Incorporated 


-16- 


Let  us  compare  this  value  with  values  which  can  be  analyzed  from 
data  obtained  by  other  investigators.  Both  Kennedy  et  al  (Ref¬ 
erence  3)  and  Schooley  (Reference  l)  have  done  experiments  on 
the  wake  of  <■  self-propelled  body  moving  in  a  fluid  with  a  verti¬ 
cal  density  gradient.  Prom  their  data  the  following  results  can 
be  computed  for  comparison. 


Kennedy  et  al 


Schooley  and  Stewart 


Slope  at  t  _  =  0.6l  ips 
col 


1.22  ips 


a  =  0.003  ft 


Slope  at  t  1  =  2.71  cm/sec 


=  5-42  cm/sec 
a  =  0.1585  ft"1 


max 


2.16 


=  4.32 


max 


3  in 
6.5  in 


=  1.1  cm 
=  4.75  cm 


Vai 


=  0.31  sec 


Vai 


=  2.26  sec 


=  0.604 


'max 


=  0.504 


y  Vag 
^ma;:  ¥ 


tooi  =  3'04 


t,nl  VaF  =  2.26 


=  1.2 


^max^o’ 


=  1.26 


HYDRONAUTICS,  Incorporated 


-17- 

Stockhausen,  Clark  and  Kennedy  deduced  the  size  of  the  mixed 
region  from  isochlor  maps  obtained  from  measurements  using  con¬ 
ductivity  probes.  With  the  occurrence  of  internal  waves,  however, 
this  method  is  a  rather  inaccurate  one.  A  plot  of  the  heights 
and  widths  of  the  envelopes  of  the  zone  of  disruption  versus 
distance  behind  the  model  was  estimated,  from  which  the  above- 
mentioned  values  could  be  computed.  In  spite  of  the  inaccurate 
way  of  determining  the  height  and  width  of  the  wake,  the  results 
are  not  too  different  from  those  in  the  present  study.  The  di¬ 
mensionless  time  of  collapse  will  be  somewhat  larger  for  the  case 
with  a  self-propelled  body  because  at  t  =  0  in  the  present  study 
the  wake  already  had  a  certain  finite  size. 

Schooley  and  Stewart  used  a  much  larger  density  gradient 
than  in  our  studies.  In  relation  to  the  large  a-value,  the  value 
y 

for  ^ax  seems  rather  large,  but  is  due  to  the  high  turbulence 
yo 

velocities.  The  dimensionless  time  of  collapse  is  larger  again 
than  in  the  present  study,  as  has  been  explained  above.  The 

value  of  1.26  for  - -  agrees  well  with  the  value  1.25 

in  the  present  study.  In  both  studies  the  value  for 

ymax  V^S 

is  slightly  higher  (by  about  50  percent)  than  in  the  present  case. 
The  reason  may  be  sought  in  the  actual  density  distribution  in¬ 
side  the  wake.  Figure  13  shows  the  difference  between  the  theo¬ 
retically  assumed  density  distribution  and  a  possible  actual 
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density  distribution.  Of  course,  the  actual  density  distribution 
in  a  wake  behind  a  self-propelled  body  may  be  somewhat  different 
from  that  resulting  from  the  use  of  a  spiral  paddle.  The  degree 
of  mixing  in  the  wake  is  a  problem  which  is  very  important  when 
one  starts  to  compare  one  experimental  study  with  another.  The 
effect  of  the  degree  of  mixing  should  be  studied  further  in  the 
future. 

Townsend,  in  Reference  7,  states  a  criterion  for  the  exis¬ 
tence  of  turbulent  motion  in  the  presence  of  a  sharp  density 
interface.  Townsend  checked  his  considerations  with  a  series  of 
experiments  in  which  a  liquid  Jet  was  injected  horizontally  along 
the  interface  between  a  dense  solution  of  salt  in  water.  From  a 
cine  record,  estimates  could  be  made  of  the  velocity  and  width 
of  the  Jet.  It  was  found  that  entrainment  of  fluid  by  the  Jet 
almost  ceased  when  the  Richardson  number  (which  increases  with 
distance  from  the  nozzle)  exceeded  0.3.  Beyond  this  point,  the 
velocity  and  cross-section  of  the  Jet  remained  nearly  constant, 
indicating  very  little  entrainment  and  presumably  very  little 
turbulent  motion.  The  Richardson  number  criterion,  then,  de¬ 
termines  the  conditions  under  which  the  stability  of  the  water 
reaches  a  sufficient  value  so  that  turbulence  cannot  act  against 
it  and  vertical  motion  ceases  entirely  (see  also  Reference  8). 

Townsend  computed  the  critical  Richardson  number  in  his 
case  as 


=  0.064  g  [(pi-p2)/p]  ( D/TJmE ) 


[12] 
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where 


D  was  the  observed  width, and 

U  the  mean  velocity  of  advance  .  i 
m 

When  applying  this  to  the  present  study,  this  Richardson  number 
becomes 


R^  *  0.064  g 
c 


by  ymax  Po  2ymax  u'2 


[13] 


Substituting  Equation  [11]  gives 


k,  =  0.128 


O.36 


[14] 


Prom  data  given  by  Townsend  (Reference  5)  the  ratio  of  u'2  and 


U  can  be  estimated  as 
m 


0.22  and 


0.45 


max 


max 


[15] 


The  above  ratios  combined  and  squared  results  in 


=  0.24 


[16] 
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The  critical  Richardson  number  for  the  present  study  is  then  com¬ 
puted  as 

=  0.128  x  7-7  x  0.24  =  0.24  [17] 

c 

Townsend  (Reference  7)  found  experimentally  the  value  0.3  as  the 
critical  Richardson  number.  This  agreement  is  remarkably  close. 

Horizontal  Spreading 

With  respect  to  the  horizontal  extent  of  the  wake  (see  Fig¬ 
ure  12),  Kennedy  stated  that  the  horizontal  extent  appears  to 
increase  almost  linearly  with  time.  From  Figure  12  it  is  clear 
in  our  case,  too,  that  initially  the  horizontal  extent  of  the 
wake  grows  linearly  with  time.  However,  as  time  goes  on,  the 
curve  becomes  convex.  During  the  time  of  collapse,  the  curvature 
changes  and  becomes  concave  for  a  short  time  interval,  then  changes 
again  into  a  convex  curve.  Finally,  the  relationship  becomes  more 
or  less  linear;  again  afterwards,  of  course,  the  curve  becomes 
once  more  convex  since  due  to  inadequate  mixing,  the  wake  never 
extended  horizontally  all  the  way  to  the  sidewalls  of  the  tank. 

All  these  changes  in  curvature  are  demonstrated  by  one  example 
given  in  Figure  14.  In  Figure  12  the  initial  changes  can  be  ob¬ 
served  more  accurately  for  several  other  runs.  The  above-men¬ 
tioned  changes  in  curvature  can  only  be  observed  if  many  data 
points  are  taken.  If  there  are  only  several  data  points  or  one 
neglects  small  changes  in  curvatures,  then  the  curve  can  be  ap¬ 
proximated  by: 


HYDRONAUTI CS,  Incorporated 


-21- 

1.  A  straightllne  section  through  the  first -measured 
points, 

2.  A  convex  part  of  the  curve,  changing  into 

3.  Another  straightline  section,  and 

4.  Once  more  a  convex  part. 

This  configuration  can  be  seen  from  Figure  14.  Moreover  the  curve 
of  horizontal  width  of  the  wake  versus  time  obtained  by  Schooley 
and  Stewart  has  these  same  general  features. 

Collapsing  Stages 

With  data  about  the  horizontal  extent  of  the  wake  available, 
a  comparison  can  be  made  with  the  experiments  of  Wu  (Reference  3). 
With  no  turbulence  present  in  his  experiments  he  defined  three 
processes  of  collapse:  initial,  principal,  and  final  stages.  He 
derived  empirical  formulae  to  describe  the  wake  collapse  of  the 
first  two  stages.  The  initial  collapse  stage  he  expressed  as 
follows: 


x-x 

0 

-  -  c 

X 

o 


[18] 


where  xq  is  the  horizontal  extent  at  the  moment  of  collapse. 


and  t'  the  time  after  collapse.  The  constants 


"c"  and 


x-x 

be  found  when  - - 


is  plotted  versus  t 


*£  g 
^  S' 


it it 

n  can 


For  the 
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present  study,  this  is  carried  out  for  three  experiments  of 
Phase  I,  one  of  Phase  II  and  one  of  Phase  III  (see  Figure  15 ). 

The  average  slope  of  the  straightline  segments  is  0.9^  and  the 
average  dimensionless  time  reading  for  their  intersections  on  the 


line 


x-x 
_ o 

X 


0.1 


\  o  / 

be  determined  to  be: 


is  0.6. 


The  values  for  Mc"  and  "n"  can  now 


Present  study 

c  =  0.16 
n  =  0.9^ 


Wu 1 s  experiments 

c  =  0.29 
n  =  1.08 


In  order  to  find  the  relationship  for  the  principal  collapse 


—  is  plotted  versus 
X 

o 


t *  y  fy  s‘  The  straightline  seg¬ 


ments  shown  in  Figure  16  have  an  average  slope  of  0.34,  while 


the  average  of  their  intersection  with  the  line  —  =1  is  at 

o 


g  =  0.95*  The  expression  for  the  principal  stage  of 
collapse  can  then  be  given  as 


=  c 


[19] 
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where  the  values  for  the  constants  are  c  =  1.02  and  n  =  0. 
while  Wu  found  c  =  1.03  and  n  =  0.55. 

The  reason  that  the  "c"  and  "n"  values  in  the  present  study 
are  different  from  the  values  Wu  found  is  undoubtedly  that  in  his 
experiments,  the  gravitational  effect  was  made  to  be  predominant 
throughout  the  motion.  In  the  present  study,  however,  both  tur¬ 
bulence  and  gravity  play  an  important  role.  Another  reason  is 
that  in  Wu’s  experiments  the  fluid  inside  the  wake  was  more  com¬ 
pletely  mixed  and  had  a  density  very  close  to  pQ.  As  has  been 
mentioned  before, the  density  distribution  inside  the  wake  looks 
approximately  like  the  distribution  given  in  Figure  13 . 

The  above-mentioned  reasons  also  help  explain  the  difference 
in  area  measurements  of  the  wake  in  both  studies.  Wu  found,  in 
the  absence  of  turbulence,  no  change  in  the  wake  volume  during 
collapse,  and  thus  concluded  that  the  mixed  fluid  inside  the 
wake  maintains  its  density  during  the  wake  collapse.  As  a  com¬ 
parison,  the  areas  for  several  tests  in  the  present  study  have 
been  plotted  in  a  convenient  way  in  Figure  17.  Before  collapse 
the  area  increases  linearly  with  time.  Then  the  curve  becomes 
convex  and  undergoes  several  undulations,  presumably  due  to 
internal  waves.  The  area  of  the  wake  still  seems  to  increase 
considerably  after  the  wake  has  collapsed.  Besides  these  two 


reasons : 
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1.  The  presence  of  turbulence,  and 

2.  The  density  of  fluid  inside  the  wake  is  not  con¬ 
stant,  there  is  a  third  important  reason  (experimental)  for  the 
apparent  area  increase: 

3.  As  time  increases,  the  wake  created  by  the  agitation 
of  the  paddle  becomes  less  two-dimensional  so  that  the  actual 
rate  of  growth  in  volume  may  be  less  than  the  apparent  rate  as 
observed  through  the  front  wall. 

CONCLUSIONS 

The  growth  of  a  wake  in  a  density-stratified  medium  can  be 
studied  in  detail  by  means  of  a  spiral  paddle  agitated  by  a 
pendulum-type  arrangement.  In  the  vertical  direction  the  rate  of 
growth  of  the  wake  is  Initially  constant,  but  finally  ceases, 
whereupon,  the  wake  collapses  vertically.  This  initial  rate  of 
growth,  the  maximum  vertical  extent  of  the  wake,  and  the  time  of 
collapse  have  been  related  to  the  Vaisala  frequency.  Independent 
of  the  size  of  the  spiral  paddle  ard  the  type  of  pendulum  arrange¬ 
ment,  two  empirical  constants  have  been  determined: 


bi 


g 


=  1.25 
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and 


t 


col 


g  =  1.25 


Besides  accurate  measurements  of  the  vertical  extent  of  the 
wake  a  precise  description  is  also  given  of  the  horizontal  spread- 
ing  of  the  wake.  These  measurements  were  used  to  estimate  the 
turbulence  intensity  at  the  time  of  collapse  in  the  presence  of 
a  density  gradient.  It  was  found  that  the  ratio  of  the  kinetic 
energy  and  the  gain  of  potential  energy  at  the  moment  of  collapse 
is  a  constant  which  may  be  expressed  as  follows: 


=  O.36 


There  is  a  need  for  a  parameter  to  take  into  account  the  de¬ 
gree  of  mixing  inside  the  wake.  This  is  especially  important 
when  comparing  different  experimentally  obtained  data.  After 
collapse,  the  horizontal  spreading  was  analyzed  for  several  runs. 
The  occurrence  of  three  stages  of  collapse  as  observed  by  Wu,  was 
checked  in  the  present  study  where  gravitational  forces  as  well 
as  turbulence  play  an  important  role.  All  three  states  were 
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clearly  noticeable.  The  constants  in  the  empirical  formulae 
•describing  the  initial  and  principal  stages  of  collapse  were  found 
to  be  somewhat  different  from  those  found  by  Wu,  due  to  incomplete 
mixing  of  the  fluid  inside  the  wake  and  the  presence  of  turbulence. 

The  volume  of  the  wake  after  collapse  seemed  to  increase 
further,  possibly  as  a  consequence  of  the  presence  of  turbulence. 
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TABLE  1 

Experimental  Results 


inch  inch  sec  ips 


— m 

0.197 

0.272 

0.211 

0.30 

0.274 

0.272 

0.355 

0.403 

0.355 

0.378 

0.456 

ragi 

0.584 

C.616 

EH 

Phase  I 


Phase  II 


0.000 

0.002 

0.00273 

0.00273 

0.00542 

0.00537 

0.C1074 

0.0210 

0.0203 


0.0000 

0.00545 

0.00547 

0.0109 

0.0217 


1.9 

1.825 

1.56 

1.71 

1.50 

1.54 

1.61 

1.60 
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FIGURE  1  -  DETAIL  OF  SALTWATER  DIFFUSER 


FIGURE  2  -  SPIRAL  PADDLES,  4  IN  AND  2  IN  DIAMETERS 
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t  -  0.25  SEC 


t  =  1 .50  SEC 


t  =  3.5  SEC  t  =  7  SEC 

FIGURE  3  -  SAMPLE  PICTURES  OF  SPREADING  AND  COLLAPSE  OF  WAKE 
FOR  a  =  0.0054  (FRONT  VIEW) 
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FIGURE  4  -  DEFINITION  SKETCH  OF  WAKE 


VERTICAL  ORDINATE  y,  INCHES 
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TIME  t,  SECONDS 


FIGURE  5-  GROWTH  AND  SUBSEQUENT  COLLAPSE  IN  THE  VERTICAL 
DIRECTION  -  PHASE  I 


X-X 


HYDRONAUTICS,  INCORPORATED 


VERTICAL  ORDINATE  y,  INCHES 
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FIGURE  6  -  GROWTH  AND  SUBSEQUENT  COLLAPSE  IN  THE  VERTICAL 
DIRECTION  -  PHASE  II 
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FIGURE  6a  -  DIMENSIONLESS  TIME  HISTORY  OF  VERTICAL  ORDINATE  OF 
WAKE  -  PHASE  II 
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TIME  t,  SECONDS 


FIGURE  7  -  GROWTH  AND  SUBSEQUENT  COLLAPSE  IN  THE  VERTICAL 
DIRECTION  -  PHASE  ill 
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FIGURE  7a  -  DIMENSIONLESS  TIME  HISTORY  OF  VERTICAL  ORDINATE 
OF  WAKE  -  PHASE  III 
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FIGURE  8  -  DIMENSIONLESS  REPRESENTATION  OF  DATA 
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FIGURE  9  -  RELATION  BETWEEN  MAXIMUM  VERTICAL  EXTENT  OF  WAKE 
AND  TIME  OF  COLLAPSE 
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o  2"  PADDLE,  LONG  PENDULUM 
A  2"  PADDLE,  SHORT  PENDULUM 
□  2"  PADDLE,  LONG  PENDULUM 
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FIGURE  11a  -  WAKE  SHAPE  DURING  GENERATION 


FIGURE  lib  -  WAKE  SHAPE  DURING  GROWTH 


FIGURE  11c  -  WAKE  SHAPE  AT  TIME  OF  COLLAPSE 


FIGURE  lid  -  WAKE  SHAPE  AFTER  COLLAPSE 


HYDRON  AUTICS,  INCORPORATED 


1.0 


0  1  2  3  4  5 


TIME  t,  SECONDS 


FIGURE  12  -  EFFECT  OF  COLLAPSE  ON  RATE  OF  SPREADING  IN 
THE  HORIZONTAL  DIRECTION 
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FIGURE  13  -  COMPARISON  OF  THEORETICAL  AND  ACTUAL  DENSITY  DISTRIBUTION 
INSIDE  THE  WAKE 


HYDRONAUTICS,  INCORPORATED 


FIGURE  14  -  HORIZONTAL  SPREADING  FOR  a  *  0.002,  PHASE  I 
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FIGURE  16  -  PRINCIPAL  STAGE  OF  WAKE  COLLAPSE 
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FIGURE  17  -  INCREASING  AREAS  OF  WAKE  WITH  TIME 
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